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ABSTRACT 



After having reported the detection of X-rays emitted by the peculiar system 
HD 5980, we assess here the origin of this high-energy emission from additional 
X-ray observations obtained with XMM-iVewion. This research provides the first 
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detection of apparently periodic X-ray emission from hot gas produced by the 
collision of winds in an evolved massive binary outside the Milky Way. It also 
provides the first X-ray monitoring of a Luminous Blue Variable only years after 
its eruption and shows that the dominant source of the X-rays is not associated 
with the ejecta. 

Subject headings: stars: individual (HD 5980) - stars: Wolf-Rayet - stars: winds 
- Magellanic Clouds - X-rays: binaries 



1. Introduction 



The most massive stars, which are also the most luminous (and for most of their lives 
the hottest) non-degenerate stellar objects, have a large impact on their host galaxies. For 
example, throughout their evolution, these stars convert simple into more complex elements, 
and distribute them into space by a complex process of mass loss. T hey are thus largely 
responsible for the chemical enrichment of the Universe (IMassey 1120031 ). Their intense, hot 
photon emission is also able to ionize the surrounding medium (leading to the formation of 
Hll regions) and to drive a continuous ejection of matter in the form of a stellar wind. Finally, 
massive stars inject a large amount of mechanical energy into their host galaxies, through 
e.g. supersonic winds, large transient erupt ions (like thos e of Luminous Blue Variables, e.g. 
i] Carinae), or GRB/supernova explosions (IMassey 1 120031 ) . 



In massive early-type binaries, the supersonic outflow from one star collides with that of 
its companion. This collisi on provokes a strong heating of the shocked wind, leading to the 
emission of X-rays (see e.g. IStevens et al. Ill992l ). Hints of this high-energy phenomenon were 
already found t wo decades ago , when the first X-ray observator ies detected overluminosities 
in hot binaries (jPollock 1 119871 ; IChlebowski fc Garmany I Il99ll ). However, since additional 
X-ray emission can have several origins, not restricted to colliding winds (CW), indisputable 
observational evidence of CW X-ray emission had to await the advent of a new generation of 
sensitive X-ray facilities (ROSAT, Chandra, XMM-Newton) . Indeed, only detailed spectro- 
scopic investigations and careful monitoring in the X-rays can bring to light the actual proper- 
ties of the hot gas produced by the wind-wind collision. Most notably, phase-locked variations 
in the X-ray domain are produced by varying separation in eccentric binaries (that changes 
the intrinsic strength of the collision) or varying line-of-sight opacity as the stars revolve 



1 Based on observations collected with XMM-#ewton, an ESA Science Mission with instruments and 
contributions directly funded by ESA Member States and the USA (NASA). 
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around each other: see e.g. HD 152248 (0 7.5III+07III, P=6 d, e=0.13: ISana et al. 1120041 ) 



=0.33; 



Gamen et al. 


2006 


; Schild et al. 


2004) 



HD 93 403 (05.5IH+07V, P=15d, e=0.23: iRauw et a 

e 
e 



20021 ). WR25 (WN6+0, P=208d, 
Pollock fc Corcoran 1120061 ). and 7 2 Vel (WC7+07.5III, P=79d, 



Up to now, all studied X-ray colliding wind (XCW) systems 
belong to our Galaxy - a few XCW candi dates have been proposed in 3 Doradus, but solely 
on the basis of an X-ray overluminosity (IPortegies Zwart et al. 1 120021 ) . which is insufficient 
to ascertain the true nature of the emission. Studying the CW phenomenon in other galax- 
ies can provide an important probe of the mass-loss process in different environments with 
different metallicities. 

HD5980, the most peculiar massive star in the Small Magellanic Cloud (SMC), lies on 
the periphery of the large cluster NGC346 associated with the giant Hll region N66. It is a 
multiple system whose main component, star A, underw ent two LBV-like eruptio ns in 1993- 
1994, increasing its brightness by up to 3 magnitudes (I Jones fc Sterken II1997I ). Together 
with star B, believed to be an early Wolf-Rayet star of the nitrogen sequence, it forms a 
close eclipsin g binary system whose perio d is 19.3d, eccentricity is 0.3 and inclination very 
close to 90° (jSterken &: Breysacher 1 119971 ) . Using the ephemeris of the latter authors, star 
A (resp. B) eclipses its companion at 0=0 (resp. 0.36), whereas periastron and apastron 
occur at phases 0.09 and 0.59, respectively. A third stellar object, star C (probably an early 
O-type star), contaminates the light of the system, but it is still unclear if this star is just 
a line-of-sight coincidence or an object gravitationally bound to the close AB pair. The 
presence of an additional component (e.g. a close orbiting neutron star) has been proposed, 
based on a 7-hour periodicity seen in spectral an d photometric variations and stochastic 



polarimetry changes, but is still debated (see e.g. IVillar-Sbam et al. II2003I ). iMoffat et al. 



(119981 ) infer the presence of a wind- wind collision region in the AB system on t he ba sis of 
strong emission - line y ariability. From UV observations, iKoenigsberger et al. I (120001 ) and 
Koenigsberger I (120041 ) conclude that the orientation of the shock cone is such that it wraps 
around star B, consistent with the notion that star A possesses the more powerful of the two 
winds. 

The de t ection of X-ray emission from HD5980 was reported for the first time by 



Naze et al. 


(2002) 


(Naze et al. 


2002. 



(2002) . Sho rt-term and long-term variations in X-ray brightness were detected 



20041 ) 



However, these observations were not sufficient to determine the 
nature of these variations (e.g. variations driven by changes in the shocked ejecta from 
the 1994 eruption, or phase-dependent changes due to wind-wind collisions in the AB bi- 
nary), and the origin of the X-ray emission therefore remained uncertain. To resolve this, 
we acquired additional monitoring of HD5980 with XMM-iVewion. 



This letter is organized as follows. Section 2 describes the observations while Sections 3 
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and 4 present the results of the dedicated XMM -iVewfom campaign and their interpretation, 
respectively. Section 5 gives our conclusions. 



2. Observations and Data Reduction 

During the past few years, HD 5980 was observed once with Chandra during 100 ks and 
five timesj with XMM-iVewfon for approximately 20 ks each time (see Tabled]). For the 
latter observations, the three European Photon Imaging Cameras (EPICs) were operated in 
the standard, full-frame mode (except for the first pn dataset) and a medium filter was used 
to reject optical light. The observations sampled crucial phases of the 19.3d orbit (Fig. [TJ; 
moreover, the last three observations were obtained during the same 19.3d orbit of HD 5980. 
To ensure a coherent reduction, we used the Science Analysis System (SAS) software, ver- 
sion 7.0, to re-process all X-ray data. After the pi peline chains, filters recommended by the 



SAS team were applied (see e.g. iNaze et al. II2006I ). For consistency, no additional temporal 
filtering was done. 

To each XMM-JVewfon dataset, we applied the SAS source detection algorithm (ede- 
tect_chain) in a region of 150" radius around HD 5980, in order to derive the best value of the 
centroid of the X-ray emission. Since HD 5980 a ppears surrounded by a soft, X-ray bright 



supernova remnant (SNR, see INaze et al. 1 120021 ) . the detection was restricted to the 1.5 



10. keV range, to minimize the contamination from this soft extended source. As a check, 
we compared the derived value of the count rate with the output of another algorithm (SAS 
task eregionanalyse) that simply calculates the number of counts in the region considered^] 
and corrects them only for the Encircled Energy Fraction (about 0.8 in our case). Within 
the error bars, both methods agree with each other so we only present the results from the 
latter (see Table [T]). 



3. Results 

It must first be noted that two XMM.-Newton datasets were taken, intentionally, at 
the same phase (<f) ~ 0.36). Although these two observations were obtained during XMM- 



2 A sixth observation taken during Rev. 0970 was strongly affected by a flare, rendering the data unusable. 

3 We used a 25" region centered on HD 5980, together with a background region of the same area but 
offset from HD 5980 by 22s east in RA (or -2040 px in X) and 82" south in DEC (or -1640 px in Y). Using 
an annular background region yields the same results, though noisier. 
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Newton Revs. 0157 (2000 Oct.) and 1094 (2005 Nov.), i.e. separated by 5 years or 97 
orbital revolutions of the AB system, they present very similar values for the count rate of 
HD 5980 (Fig. [2] and Tabled]). This means that post-eruption, epoch-dependent variations 
are minimal. Indeed, if a significant part of the X-ray emission came from the collision of the 
fast wind with the slower ejecta associated with the 1994 eruption, a monotonic decrease in 
the X-ray brightness and a decline in the X-ray temperature with time would be expected as 
this interaction zone expands and cools. The signature of this slow/fast wind interaction zone 
was detected as discrete Siiv absorption components at —680 km s^ 1 i n 1999, only one year 



before the fi rst X-ray observation, accelerating to —710 km s in 2000 (IKoenigsberger et al 



2000l . l200ll ). and disappearing by 2002. Our XMM-jVewion observations however indicate 



that any observable X-ray emission in the 0.5-10 keV band from this interaction is minimal. 

Second, coherent variations with phase are now clearly detected, with X-ray emission 
steadily increasing towards <fi = 0.36 (Figs. [2]and[3]). At that time, the emission also appears 
slightly harder, though we cannot exclude a constant value of the HR at the 2a level. This 
bright phase corresponds to the eclipse of star A by its companion, i.e. when the opening 
of the shock cone, slightly concave with respect to star B, is directed towards the observer 
(Fig. EJ. 

Although we cannot directly compare the Chandra and XM.M.-Newton total count rates 
(because of background contamination, different spatial resolution, and cross- calibration 
problems), it must be noted that the variations measured by Chandra confirm the trend 
seen by XMM-iVewion. During the 100 ks of the Chandra ob servation, the coun t rate clearly 



increased from 2 to 4xl0 _3 ctss _1 in the 0.3-10. keV band (INaze et al. 1 120021 ) or from 1 to 
2xl0 _3 ctss _1 in the 1.5-10. keV band (Fig. [2]). These data covered phases ranging from 
0.24 to 0.30. In this interval, the X.M.M.-Newton EPIC-MOS lightcurve predicts an increase 
of ~ 1 x 10 _3 ctss~ 1 in the count rate. Simulations made with PIMMqfl predict the EPIC 
MOS and Chandra count rates to be similar in this energy range for a large range of spectral 
parameters: the variations observed by the two X-ray facilities are thus fully compatible. 

Finally, we investigated the X-ray spectrum of HD 5980. Of course, the faintness of 
this distant source and the contamination at low energies by the superposed SNR prohibit 
any detailed study. However, since the spectral propertie s of the SNR conta mination can be 



determined thanks to the high-resolution Chandra data of INaze et al. I (120021 ). a first, general 
spectral evaluation can be made: our data only reveal a clear, usable excess at high energies 
for Revs. 0157, 1093 and 1094 (i.e. when HD 5980 is the brightest). This high-energy 
emission was fitted by an absorbed optically-thin plasma model (mekal, in Xspec vll.2.0): 



4 ht tp : / /heasarc . gsfc . nasa.gov/Tools / w3pimms . html 
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compared to the spectral properties reported by lNaze et al. I (120021 ). only changes of the flux 
level were detected; within the confidence intervals, the absorption and temperature do not 
seem to vary in a significant way. However, it must be noted that the spectra are rather 
noisy and that only large variations would be detected. 



Discussion 



The observed variations of the X-ray emission from HD5980 can be qualitatively ex- 
plained by considering the geometry of the CW region in the close A+B binary. 

The change in the X-ray flux could be related to absorption. Indeed, at phases close to 
0.36, the opening of the shock cone is directed towards the observer an d the X-ray emission 
of the CW region is thus seen through the lower density wind of star B (IKoenigsberger 112004 
and Fig. [T]). The lower absorption results in an increase of the observed X-ray luminosity. 
At other phases, the very dense, slower wind of star A strongly absorbs at least part of the 
CW emission. Such a scenario h as been proposed to explain the behaviour of the WR+O 
binary 7 2 Vel flSchild et al. 1 120041 ) . If the X-ray li ghtcurve of HD 5980 is symmetric with a 
peak at 0=0.36, the width of the lightcurve can be estimated to FWHM=0 .26 in phase, 



which corresponds to a half opening angle of the shock c one 9 of about 46° (jWillis et al 



19951 ) and a momentum ratio of 3-4 (jStevens et al. 1119921 ). suggesting that the wind of star 



A is even stronger than shown in Fig. [IJ 

However, the available observations do not yet enable us to ascertain that the maximum 
brightness really occurs at/near 0=0.36. For example, it is conceivable that the count 
rate actually continues to increase beyond this phase, e.g. to peak at apastron. A similar 
variation, proportional to the separation between the stars rather than inversely proportio nal, 



was observed for Cyg OB2 #8A (06I+05.5III, P=22d, e=0.24; iDe Becker et al 



20061 ). In 



that case, the variation is explained by a sudden decrease of the wind velocity at periastron 
because of a strong radiative inhibition/braking or as a result of perturbations associated 
with tidal bulge interactions. For HD 5980, such changes of the wind velocity would mean 
that the X-ray emission should not only be stronger, but also harder towards apastron, as 
suggested by the available data. 

To distinguish between these two scenarios, additional X-ray observations are necessary 
in order to complete the lightcurve near apastron. This would enable us to notably determine 
the position of the brightness peak (apastron or eclipse of A by B ?) and the exact shape 
of this increase. A detailed spectral analysis is also needed, but it requires better spatial 
resolution than that of XMM -iVewfon in order to cleanly disentangle HD5980 from the 
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superposed SNR, therefore permitting access to the lower energy range which is more affected 
by the absorption variations. 

In addition, hydrodynamical modelling of the fast wind/slow wind interaction region 
and the CW area should be undertaken to quantify their respective strength and expected 
modulation. In this context, the lack of strong secular variations over 5 years is rather 
puzzling, since, during this interval , the wind density of star A changed, as witnessed in 
the UV and optical emission lines (IKoenigsberger 1 12004) . Howeve r, the X-ray luminosity 
(in the radiative limit) goes as Lx ocMxv^ ( iStevens et al. Ill992l ). and since the eruption 
the terminal v elocity of the w i nd of star A has increased from about 600 km s _1 to about 
2000 km s _1 ( IKoenigsberger 1 12004 ). which offsets the order-of- magnitude decline in the 
mass loss rate over that interval. On the other hand, the relative constancy of the X-ray 
brightness at <fi = 0.36 despite a significant decline in the mass loss from star A might indicate 
that star A's wind has mainly changed in the direction perpendicular to the orbital plane 
than it has in the orbi tal plane, as might be ex pected if the mass loss from star A is not 
spherically symmetric (jVillar-Sbaffi et al. 1 120031 ) . Note however that the observed P-Cygni 
UV absorptio n lines indicate that a t least some of star A's wind variation occurred near the 
orbital plane (IKoenigsberger 112004 , and references therein). 



5. Summary and Conclusions 

An X.M.M.-Newton monitoring campaign of the peculiar massive binary HD 5980 has 
unveiled the variations of its X-ray emission. As two observations taken 5 years apart (a 
separation of more than 90 binary orbits!) present similar count rates, the main source 
of X-rays cannot be the fast wind/slow wind interaction following the LBV-like eruption 
of star A, since its emission is expected to monotonically and rapidly decrease with time. 
Because individual Wolf-Rayet stars and LBVs are only weak X-ray sources, the high-energy 
radiation must be associated with the wind-wind collision in the close A+B pair, a fact 
further supported by the detection of phase-dependent changes. This is the first time that 
the presence of X-ray emitting gas produced by the collision of winds in a binary has been 
confirmed outside our Galaxy. 

The X-ray emission of HD 5980 appears modulated with phase: a clear increase is 
observed towards the time of the eclipse of star A by its companion. This can be due either 
to the lower absorption inside the shock cone (a situation reminiscent to that of 7 2 Vel) or to 
a lower wind velocity at periastron (similar to the behaviour of Cyg OB2 #8A). Additional 
data and hydrodynamical modelling are now needed to further distinguish between these 
two scenarios. 
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Table 1: Results from the XM.M.- Newton observing campaign. Count rates were evalu- 
ated by the simple eregionalyse task and are given in the 1.5-10. keV band (in units of 
10" 3 ctss _1 ). Hardness Ratios (HR) are defined as (H - M)/(H + M) with M and H the 
count rates in the 1-2 and 2 -10 keV bands, respectively. Phases refer to the ephemeris of 



Sterken &: Breysacher I (119971 ). while dates are Julian Dates minus 2 450 000d. 



Rev. 


Date 


Duration 


<t> 


Count rate MOS1 


Count rate MOS2 


Count rate pn 


HR (MOS1) 


0157 


1835.245 


0.245d 


0.37 


8.26±0.85 


6.95±0.79 


14.92±1.47 


-0.46 ±0.06 


0357 


2234.643 


0.322d 


0.10 


3.00±0.54 


2.54±0.60 


6.41±1.16 


-0.77 ±0.08 


1093 


3701.863 


0.218d 


0.26 


4.64±0.67 


5.16±0.75 


17.32±1.44 


-0.60 ±0.07 


1094 


3703.807 


0.207d 


0.36 


6.13±0.80 


8.14±0.88 


17.41±1.45 


-0.56 ±0.07 


1100 


3716.125 


0.208d 


1.00 


1.70±0.73 


2.60±0.69 


8.38±1.29 


-0.80 ±0.10 
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Distance (R ) 



Fig. 1. — Representation of geometry of HD 5980's wind-wind collision region at the 4 orbital 
phases for which XM.M-Newton data are avail able. The shape of the shock surfaces was 
computed under the thin shock approximation (ICanto et al. Ill996l ) and assuming that Ma 
= lxl(T 4 M yr" 1 , M B = 2xl(T 5 M yr" 1 , v oo (A)=2000 km s -1 , and v oo (B)=2600 km 



s 1 . The broken circle arou nd star A represents the extent of the wind accelerating region 



(IKoenigsberger et al. II2006I ). 




Fig. 2. — Bottom: Evolution of the count rates and hardness ratios with phase from XMM- 
Newton observations (see definitions in Tabled]). Top: Count rate of HD 5980 measured 
during the 100 ks Chandra observation, in the same energy band. 
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Fig. 3. — Combined MOS images of the area surrounding HD 5980 in the 0.5-10. keV energy 
range. The data have been binned to get pixels of 2.5" and then smoothed with a gaussian 
of a— 3px. 



